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The semiclassical orbitals of a relativistic electron on a ro- 
tating sphere threaded by an intense magnetic dipole field 
are examined. Several physically distinct regimes emerge, 
depending on the relative sizes of the mass, total energy, 
canonical azimuthal angular momentum, and magnetic field 
strength. Magnetic flux enclosed by orbits is quantized very 
close to the poles, suggesting a quantum Hall-like state. Ap- 
plication of this system to neutron star surfaces is outlined. 
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1. INTRODUCTION 



Neutron star pulsars are known to have surface mag- 
netic fields of up to 10 12 ~ 13 Gauss |j]-[§- The magnetic 
field lines dragged in by neutron star collapse are pre- 
sumably squeezed by magnetic flux conservation in the 
contracting plasma. However, that assumption does not 
answer the question of how the magnetic field sustains 
itself at later times. The crust, and at least part of the 
interior of a neutron star, have electrons and protons. 
Their internal macroscopic currents might also affect the 
field, if not actually create it, but these internal currents 
are in turn strongly affected by the intense field. We 
assume that such currents are essentially electronic. If 
the density of electrons is low and/or the magnetic field 
large, a phenomenon akin to the quantum Hall effect ||] 
should be expected. 

We select a simplified system retaining some features 
of a realistic neutron star. This system is chosen to iso- 
late quantum Hall-like surface states: a charged particle 
constrained to move on the surface of a sphere of radius 
R threaded by an intense magnetic dipole. The treat- 
ment is restricted to surface states as their distinctive 
properties may someday be observable. The sphere ro- 
tates with angular velocity ft, not necessarily parallel to 
the magnetic axis. Our main discussion assumes that the 
two axes are parallel in order to estimate the effect of ro- 
tation. This is followed by the general case of tilted axes, 
which is not greatly different for a moderate rotation rate. 
The problem can be treated relativistically by use of the 
appropriate metric for the rotating frame. The surface 



orbits serve as a prelude to the three-dimensional tilted, 
rotating system, which we treat in a separate paper [pj. 
Neglecting rotation, the distinct physical regimes can 
be characterized by comparison of three dimensionless 
parameters: 

A) = \q\B Q R/(2mc 2 ) , 
e = E/mc 2 = /9o77 , 
; = 2cP 4 ,/{qB Q R 2 ) 

which are, respectively, the magnetic field strength Bq 
in rescaled units at the magnetic poles; particle energy 
E in units of the rest mass; and particle canonical az- 
imuthal angular momentum P^, in rescaled units about 
the rotation axis. For electrons on a neutron star with 
the strongest measured fields, (3q ~ (0.1)i?/(2Ac) ~ 
2.5xl0 15 , where Ac is the electron Compton wavelength 
= 2irh/(m e c). The relevant energy scale is then set not 
by mc 2 , but by r\ — e/ (3q. 

The regime of r\ Jj 1 is the ultrahigh energy case, E ~ 
10 21 eV or higher, depending on Bq and R but not on m. 
The magnetic field /?o sets the scale for P$ : when \l^\ ~ 
1, the ultrahigh P^ case, the effect of P<f, can overcome the 
inhibiting effect of the field. When r\ and 1$ are small, the 
charged particle has no allowed region of motion on the 
spherical surface except very close to the rotational and 
magnetic poles. A large l^ allows narrow regions away 
from the poles, but only n <; 1 allows charged particle 
motion over a substantial portion of the sphere. 

The rotational angular velocity is rescaled to u> — 
flR/c, which is ^ 0.1 in realistic cases [0-||. We take 
u> = 0.1 throughout as illustrative, being an order-of- 
magnitude upper limit on pulsar rotation Q. 

The simplest treatments, in the limit of infinite field 
strength, of charged particles trapped on spheres by in- 
tense magnetic fields result in particles frozen in place 
in the crust |||J]]. The treatment here adds the feature 
of expanding the particle motion in inverse powers of 
the field strength. Although electrons are stripped from 
the neutron star surface by the rotation-induced electric 
field, the bulk of electrons remain in the crust to prevent 
significant charge separation, with the surface sheathed 
by a thin space charge. The surface space charge is sta- 
bilized by the Coulomb force (with the positive crystal) 



opposing the induced electric field. Only a small fraction 
of electrons are accelerated into the stellar wind [g] . The 
application of this paper is to quantum single- and many- 
body states where radiation emission is neglected. This 
is exact for charged particles in their ground states or in 
excited one-body states unable to decay by Pauli exclu- 
sion blocking in the presence of other fermions. We also 
seek a general classification of possible orbitals, based on 
the kinematic parameters e, rj, and 1$ and rotation u>. 



2. CLASSICAL KINEMATICS 

We explicitly show all factors of c and, in Sect. IV, of fr. 
Thus x^ = (xo,x) = (ci,x). The metric has dimensions 
with signature (H ). 



2.1 Generally covariant Lagrangian 

The general relativistic action of a particle of mass m 
and charge q in an external electromagnetic and gravita- 
tional held is 

S = -mc I ^Jg^ iv {x)dx^dx v + - A tl (x)dx fi (1) 



with fixed endpoints in parameter-independent form 
The path parameter is the proper time r. Then if 



x» 



dx^_ 
~dV 



(2) 



the general relativistic Lagrangian of a particle in the 
external fields A^ and g^ v is 



2.2 Two-dimensional rotating sphere 

We choose axes so that the magnetic dipole is along 
the 9 = direction. The threading dipole magnetic field 
has polar strength Bq : 



Ag = 

A r = Q 

A = {B R 2 /2) sin 2 9 , 



(8) 

(9) 

(10) 



where A^ is defined in the rotating spherical coordinates. 
The rotation axis is tilted at an angle 9$ with respect to 
the magnetic dipole in the <j> = 0, n plane (Fig. |l|). (This 
apparently strange choice of coordinates is motivated by 
the the fact that the magnetic field still dominates the 
rotational effects.) The metric in a spherical polar coor- 
dinate system (r, 9, 4>) rotating with the sphere is given 
by the line element 

ds 2 = g^dx^dx" 

= c 2 (l - u> 2 )dt 2 - r 2 (d9 2 + sin 2 6W0 2 ) - 

2crtu l psm9 dt d<fi — 2cru)e dt d9 . (11) 

The vector u is defined from the rotational angular ve- 
locity vector $7 by u> = il X r/c, with components 



cj0 = u>[cos 9q sin 9 — sin 9q cos 9 cos i 



U! S 



-uj sin 6*o sin < 



(12) 



9 2 



These are the appropriate generalizations to the case 

00 7^0 §. 



L 



mcy g^vX^x" H — A^x^ 



(3) 



The canonical momenta are given by 



p = *L 

11 d& 



meg^x" q 

== H A L 



g a px a xP 



(4) 



Using the identity 

9^ X = Si , 
Eq. (4) gives the constraint 

gT{P» - 1 A li ){P u - ?-A v ) = (mc) 2 
c c 

The equations of motion are given by 



dL dPf, 
dx^ dr 



= . 



(5) 



(6) 



(7) 



FIG. 1. The sphere of radius R, threaded by a magnetic 
dipole field M, and rotating with angular velocity fl which is 
tilted by angle #o with respect to dipole. 

The Lagrangian is expressed in terms of the proper 
time r, with x^ = dx^/dr. Eq. (7) are three equations, 
one for each momentum component, suppressing radial 
motion. Since the Lagrangian does not explicitly depend 
on i, we have 



dL 

!k= 



(13) 



which implies 



Because 

dT = 

we have also 



cLFb 
dr 



dt y/g^ {dx^/dt)(dx u /dt) 



dP6_ 
dt 







(14) 



(15) 



(16) 



Thus the energy E — Pq is conserved. The equations of 
motion for d> and 9 are non-trivial: 



dL dP 6 



d<j> 



dr 



and 



PL dPg _ 

39 dr 



(17) 



(18) 



The constraint (6) is 



c 

= (mc) , 



(19) 



including both g and g ^ terms. The contravariant 
metric components are 

g 00 = i/c 2 , g rr = -l/r 2 , 

g 89 = -(l-u> 2 e )/r 2 , 5 ^ = -(l-c 2 )/(r 2 sin 2 0), (20) 
.g 00 = -u>4,/(crsm0) , g oe = -w e /{cr) , 

where r — R for our case. Along the actual worldpath 
in spacetime, the constraint g^x^x^ = c 2 obtains; this 
condition is valid after varying the action and simplifies 
the equations of motion. 



3. CLASSICAL ANALYSIS: ZERO TILT 



The opposite sign of the square root, not shown, is also 
valid. We consider Eq. (22) in three different limits of 
physical interest, with q > 0. The q < case can be 
obtained from the q > case by reversing the sign of 1$ 
and ui. 

While the </>(£) motion is trivial, the polar motion 9{t) 
is not. The two angular motions, both periodic, decou- 
ple from one another because 1$ is conserved. Without 
radial motion, the polar motion alone is a closed one- 
dimensional system. The two periods, r^ and rg, are 
in general not equal or even commensurate: their ratio 
tq/t^ is not necessarily a rational number k^/kg, where 
(fc^, fcg) are a pair of integers with no common divisors. 
If the periods are commensurate, then the orbits can be 
arbitrarily complex, but close after a time r c i osuro = kgrg 
= k^r^; otherwise the orbits never close. We compute 
below the magnetic flux $ enclosed by a pole orbit (see 
III.C), but the flux is well-defined only if the orbit is 
closed. For very large fields, nonetheless, the variation 
of 9 is O(e/Po) and tiny in these two cases (ultrahigh 
P<p and localized pole orbits), and we define the enclosed 
flux by one complete revolution of <j> from to 2ir at the 
approximately constant polar angle # m id : 



$(0mid) = TTB R 2 [l-COs20 mid ] 

f2i:Tic\ /mcR\ „ , . 

= (--)(— )Ml-cos2tf mid ] 



(23) 



In general, this flux is macroscopically large, a product of 
two large dimensionless factors and the small elementary 
flux quantum. 

In a uniform field, the particle's cyclotron radius and 
magnetic flux enclosed are constant. The enclosed mag- 
netic flux is still conserved, as an adiabatic invariant, 
for slowly- varying fields P]Jq]. A charged particle on a 
sphere in the present configuration would, in general, 
see a rapidly varying field. But if the field is intense 
(/?0 — * oo ) and the energy small (rj <C 1), the variation 
of the particle's orbit from constant 9 is higher order in 
l//3o, and the enclosed magnetic flux is quasi-invariant 
(see also Sect. V). 



A simplified treatment of rotation assumes zero tilt, 
9 a = 0. Eq. (8-10) define the dipole field, with A^ de- 
fined in the rotating spherical coordinates. If #o = 0, the 
Lagrangian does not explicitly depend on </>, and 



dP^ 
dt 







(21) 



yielding another constant of motion, the canonical az- 
imuthal angular momentum P^, along with E. Further, 
we define 



lg = 2cP e /(qB R 2 ) 

= (l//3 o ){-l + e 2 -2^e/3 o [; 



sin 



[l/sm 2 0-w 2 ]/3^-sin 2 0] 2 } 



1/2 



(22) 



3.1 Ultrahigh energy orbits 

In the limit of very high energy, i.e., e ^> 1, Eq. (22) 
becomes 



h = vV + 2cdr l sin ° ~ t 1 ~ q2 sin e \ sin ° 



(24) 



where P^ has been neglected. That is, 77 ~ 0(1) and 
\l<p\ -C 1. In Fig. |[ we plot the right hand side versus 
cos 9 and various values of E. We find that, for values 
of energy rj < a critical value rj c (uj), there are four turn- 
ing points enclosing two distinct allowed regions for the 
particle between the poles and the equator, with two of 
the four turning points very near the poles. (These po- 
lar turning points are nonzero because of the centrifugal 



effect of P<f, and vanish as P<f, vanishes.) As we increase 
the energy, the allowed regions expand and, for r\ > r) c , 
merge into a single region occupying essentially the whole 
surface of the sphere. In the latter case, two of the turn- 
ing points merge and disappear. The turning points in 
the ultrahigh energy limit are: 



1 - 2u}i] ± y/1 - Auj-q 
2Q 2 



(25) 



Only one of these roots, sin 6L, is physical in the limit 



ft -> 0: 



if. If Q, = 0, then 77 c = 1. 



is thus O(l/Po). There are four endpoints in total, two 
for each allowed region. Fig. [| shows 9 as a function of 1$. 
The values of 1$ are restricted by the condition < |Z^| < 
1. 



FIG. 3. Ultrahigh P$ orbits: two allowed annuli with cen- 
tral angles of 6 (solid) and n — 8 (dashed) as functions of 
ld>- 



3.3 Localized pole orbits 



FIG. 2. Ultrahigh energy orbits: Pe{0) for energies r/ = 0.5 
4.0, in units of (J3oTf)(mcR). 



3.2 Ultrahigh P^ orbits 

Taking \l$\ in this case to be ~ 0(1), Eq. (22) can be 
re-expressed as 



le = (l//3b){-l + ( 2 - 2wr)Po[l$ - sin 2 6} 
[l/sin 2 9~0J 2 }p 2 [l <p -sm 2 9} 2 } 1/2 



(26) 



The ultrahigh magnetic field introduces a very large neg- 
ative quantity into the square root. We expand results 
in inverse powers of (3$. For a given P^ or 1$, the two al- 
lowed regions are very narrow in 9 and, to 0(1) in l//3o, 
given by the two distinct values of 9 : 



sin# = \/L 



(27) 



where 9 is defined over the range to tt, and i) « 1; 
that is, e has been neglected compared to 1^. Note that 
< 1$ < 1. To O(l/(3 ), the endpoints of the allowed 
regions are given by 



where 



<h. 



cos20 = cos20 + <5i, 2 //?o 



■oieTN' + ^-^tt 2 -!)] 1 / 2 
2[1/^-cj 2 ] 



(28) 



(29) 



where the subscripts (1, 2) here refer to the plus-minus of 
the square root. The angular width of the allowed region 



Eq. (22) can also be written, for P^ ^ 0, as 
le=(l/l3 ){-l + e 2 ~2Cuef3o[^ 



sm 



[l/sm 2 9~uj 2 }[3 {^ 



sm' 



9} 2 } 1/2 



(30) 



For r\ and \l^\ <§C 1, allowed regions appear only if 6 — > 
or 7r. This case can be understood by extrapolating 
the ultrahigh energy and P$ case to 77 and 1$ — > 0. The 
allowed regions are two narrow azimuthal annuli centered 
about the north and south poles (Fig. ||). 



FIG. 4. Localized pole orbits: narrow allowed annuli (dark- 
ened circles) near poles 6 = and n, shown on sphere of radius 
R. Magnetic dipole M and rotational ft axes are parallel in 
this case. 

The turning points near the poles are given by 
9 = 9 ± A9/2 
= \fl 



l±i 
2 



Wo% 



> 0) (31) 



and the same replacing 9 by it — 9. In both cases, as 
the field grows, the turning points approach zero angle 



with the rotation axis and the annular width vanishes. 
(The Q{1$) effects are too small to show in Fig. |[) The 
magnetic flux enclosed by an orbit at the rotational poles 
is given by 



$(0~O,tt) = 2irR 2 B 



-t2 



(A-kTic\ fmcR\ , 



(32) 



and the same replacing 9 by tt — 9. Note that $ oc Ia, oc 
Pcjj/q. The expression (32) for the magnetic flux does not 
have to be macroscopically large, as it consists of two 
large and one small dimensionless factors times the ele- 
mentary flux quantum. If 1$ is small enough, $ can be 
microscopic, signaling a quantum Hall-like state. 



4. SEMICLASSICAL QUANTIZATION: 
ZERO TILT 

Although we do not carry out the full quantum analy- 
sis, a semiclassical treatment brings out many of the de- 
sired features. For periodic classical orbits, semiclassical 
quantization is most easily implemented with the Wilson- 
Sommerfeld (W-S) or Bohr-Sommerfeld conditions ||. 
This procedure is valid for 9 quantization as 9 motion 
is a complete subsystem alone if 6*o = 0. The closure or 
non-closure of the classical orbits is then irrelevant. 

The condition on the validity of the W-S procedure is 
that the particle's de Broglie wavelength be much smaller 
than the length scale over which the background field 
varies. In our case, 

2nh 

IpI 

a condition that holds except where |p| — ► 0. Only in- 
finitesimally close to the 9 turning points does this condi- 
tion break down, as the three-momentum |p| is otherwise 
large and R enormous in any case. The W-S semiclassi- 
cal quantization is validated by the WKB approximation, 
with this restriction j|. 

The classical analysis already implies orbitals reminis- 
cent of a quantum lattice: quasi-free conduction bands 
(ultrahigh energies) and localized states (polar and az- 
imuthal orbitals). The azimuthal (ultrahigh Pa/) rings are 
localized in one, but not two, dimensions. They conduct 
along one direction but are trapped in the other. The 
pole orbitals are similarly localized, but are also confined 
in absolute position to be near the poles. 

The classically allowed regions are those in which the 
quantum wavefunctions are oscillatory rather than expo- 
nential. Because of the the Heisenberg uncertainty prin- 
ciple constraint, the quantization conditions, together 
with the classical spherical cyclotron relations, lead to 
a tightly constrained set of semiclassical orbitals. In a 
uniform, planar magnetic system (the Landau system), 



these orbitals would be equivalent to harmonic oscillator 
states §,0. 

We now impose W-S quantization. Since Pa, is con- 
served, the azimuthal quantization is trivial: 



Pa, = n^h 



(33) 



using the Bohr form of the W-S rule, valid for circular 
motion. For Pa, 3> h, tia, is very large and the </> motion 
is essentially classical. We assume this to be the case, 
except very near the rotational poles; in that case, the 
orbital size vanishes as Pa, — > 0, so that Pa, must be small. 



4.1 Localized pole states: flux quantization 

The 9 motion is periodic but non-circular and requires 
the alternate form of the W-S rule, 



Pg d9 = (ng + 1/2)1% 



(34) 



for turning points 9\ and 9i- The replacement n — ► n+1/2 
preserves the exact quantization for simple harmonic mo- 
tion. This integral is evaluated here in a simple rectangu- 
lar approximation (Simpson's rule) . Witheqs. (30,31,34), 
we obtain 



E ~ (mc 2 ) y/l + (ng + l/2)(hp Q /(2mcR)) 



(35) 



Fig. @ shows E/mc 2 as a function of Bq for ng = 0, 1, 2. 
The dimensions are set by the critical field B Cl the field 
strength with mc 2 of energy within a Compton cube of 

volume (2irh/mc) 3 : B c = (m 2 /Tr)Jc 5 /h 3 . Fields with 
this strength or greater introduce effects of quantum field 
theory such as vacuum polarization. The energy step size 
is controlled by the tiny ratio of Compton wavelength 
h/mc to i?, multiplied by the enormous magnetic field 
/3 - Electrons in excited states radiate until they reach 
the ground state (ng — 0). 



FIG. 5. Quantized pole states: energy E/mc 2 as 

a function of pole field strength Bo for quantum states 
ng = 0, 1, 2. For the electron, magnetic field units 
in B c = (m 2 e /TT)(c 5 /h 3 ) 1/2 = 1.3 x 10 12 Gauss and 
q 2 /(hc) = 1/137. 



The magnetic flux enclosed by a semiclassical orbit, 
being an adiabatic quasi- invariant, is quantized and is 
given by 



$(0~O,tt) ~ 2ttB R 2 6 



(36) 



4.4 Density of states 

The density of states is important in any charged 
fermion system for determining the conductivity. The 
semiclassical density of states dM/dE is M 



so that 
$(<9~0,tt) 



4ttcP 

q 



4Trhc\ / mcR 



)( — r— JPolcf, OC TI4, , (37) 



analogous to the quantum Hall effect. Any n@ depen- 
dence in $ is a correction to Eq. (37) of relative order 
O(l//3o) and arises from the breakdown of exact invari- 
ance of <&. 



4.2 Localized azimuthal rings 

In the limit of high azimuthal angular momentum 
|^| - 0(1), using eqs. (26-29,34), we obtain 



E 



{mc 2 ) yjl + (ng + l/2)(h/3o/(2mcRy/l^)) 



(38) 



The same combination of tiny ratio multiplied by very 
large /?o that occurred in Eq. (35) appears again, and 
excited states are again unstable to radiation. Note, the 
magnetic flux enclosed by a semiclassical orbit is not in- 
teresting in the quantum regime, as it is macroscopically 
large. 



4.3 Poleward conduction: quasi-continuum 

In the limit of ultrahigh energies using eqs. (24,25,34), 
we obtain 



E ~ C(ng + 1/2) (hc/R) , 



(39) 



where C rs 4, 1, or 1/2 for r\ ^, ~, or <; 1. Note that the 
energies scale as harmonic oscillator energies, where the 
frequency is set by the size of the sphere, not the mag- 
netic field. As in Sect. III. A, we have assumed P$, to be 
negligibly small. In realistic cases, hc/R <C mc 2 and E, 
so that the energy levels are very closely spaced. Such 
levels form a quasi-continuum that allows the charged 
particles to conduct almost freely along the polar direc- 
tions, subject only to external crystal resistance and ra- 
diation losses. For these ultrahigh energies, ng is so large 
that the motion is essentially classical. The energy levels 
depend only weakly on fi, for small flR/c. 

The magnetic flux <I> is not interesting in this case 
because the allowed region is a large part or all of the 
sphere. Flux quantization is not relevant because the en- 
ergy is independent of the magnetic field Bq, apart from 
the overall constant C. Only if the energy E is taken 
<S O(\q\B R/2) do B and P<p become important again; 
this is the localized pole orbit case. 



dN = 



d6 dPf, 



dP A 



(2tt^) 2 

and the differential surface area element is 
dS = R 2 sin 6 dO dcj) . 



(40) 



(41) 



Thus we have for the the number of states per unit area 
per unit azimuthal angular momentum per unit energy 



d 3 Af 



dPo 



1 



dSdP^dE dE\Air 2 h 2 R 2 sin 9 
Eq. (22) can be rewritten as 
Pe/mcR = [(E/mc 2 ) 2 - 



(42) 



(E/mc i )ai(9;P^Q)-ao(fi;P4„n)] 1 ^ , (43) 



with 



ai = 2{VL/mc 2 )[P <j} - (qB R 2 /2c)sh 



(QR/c) 
Therefore 



\P* 



= l + (l/sin 2 6»- 
- (qB R 2 /2c) sin 2 0} 2 /(mcR) 2 



^ = (mcR/Pe)[E/mc 2 



(44) 



(45) 



(46) 



•oi/2](JZ/c) 

The factor dPg/dE is shown as a function of n in Fig. O 
for j/^l <C 1, and as a function of E/mc 2 in Fig. for 1$ ~ 
0.1. For low energies, E — > mc , this function approaches 

[00 + a 2 /4] 



[l- ai -a r 1/2 {l-ai/2 



-{K/mc 2 )} , (47) 



1 — a\ — ao 

where K = E — mc 2 . In the ultrarelativistic limit, E ^S> 
mc 2 and/or rj <; 1, dPg/dE approaches a constant, R/c, 
a characteristic of two-dimensional Landau states. Note 
the large peak in dPe/dE at low e or 77, the semiclassical 
edge of the discretized quantum regime. 



FIG. 6. Density of states factor dPg/dE, in units of R/c, 
as a function of energy 77 for various values of 6 and |/^| ^C 1. 



FIG. 7. Density of states factor dPg/dE, in units of R/c, 
as a function of energy E/mc 2 for 1$ — 0.1 



In the quantum limit, 
(W = dng(- 



dP* 



2irh 



(48) 



treating P^ as continuous. Equation (42) becomes: 



d 3 M 



_ dn e f 



1 



dCdP^dE dE \2irhR 2 sint 



(49) 



where dC = Rsm9 d<j>, a unit of azimuthal arc 
length. The function dng/dE can be obtained from 
eqs. (34, 35, or 38). For the localized pole states of 
Eq. (35), take E 2 = (mc 2 ) 2 + (n g + l/2)ef, where e\ 
= e?(gS o ,i?,P ).Then 



dng 
~dE 



dn e _ IE 
~dK ~1T 



2(mc 2 + K) 



(50) 



for E 2 > (mc 2 ) 2 + e\/2 



Related to the density of states is the degeneracy fac- 
tor for a given ng level. Semiclassically, this degeneracy 
arises from shifting an orbital center and is identical to 
the degeneracy of planar Landau states g . The number 
of degenerate states available per dS per dP^ at state ng 
is 



d 2 N 
dS dP a 



(ng) 



( ei /4^n 2 ci?^) y/ng + 1/2 . (51) 



5. NON-ZERO TILT: APPROXIMATE ORBITALS 
5.1 Adiabatic quasi-invariance 

As in Sects. II and III, the energy E = Pq is exactly 
conserved, but P^ or 1$ no longer is, unless u)sm.9o is 
zero. Pg or lg was never conserved, but in with zero tilt, 
its motion was exactly separable with 1$ conserved. That 
separability is also lost if u> sin 9q ^ 0. 

When 6*o 7^ 0, the non-conservation of 1$ is apparent 
from the equation of motion (10): 



dT 



sgn(g) 
cRf3 



1 §_ur 

2 d<b 



!%£(cri>i8m9) + %£(cr0i) 



(52) 



using the general metric for #0 7^ 0. The non-conservation 
of l^ is 0(1/ Po), as well as requiring <Dsin# ¥= 0- I n the 
limit of ultra-intense magnetic field (3q, the motion is not 
qualitatively different from the zero-tilt case. The orbits 
are more involved and possibly chaotic, but so focused by 



the intense field that, to zeroth order in (1//3q), the re- 
sults of Sects. II and III with 9q = are still qualitatively 
valid. 

The approximate conservation of 1$ validates a pertur- 
bative treatment of the classical orbits. In Sect. II, where 
1$ is exactly conserved, the action variable 



Jih = 



Pi, 



(53) 



is an adiabatic invariant MKU ■ This allows for the semi- 
classical quantization of P^ as a trivial step in that case 
as well as the exactly separable treatment of Pg . We thus 
expect any effects of the non-zero tilt to be suppressed 
by powers of to and 1//3q. 



5.2 Classical motion 

Since the <j) motion to zeroth order in l//?o is periodic, 
global expressions involving the (f> coordinate (such as E, 
but not the orbit </>(t), 9(t)) can be averaged over the full 
circle <f> 6 [0, 2tt] with 1$ treated as constant, in order to 
determine the first-order corrections O]. The 0(1/ (3q) 
correction to dl<p/dr = is obtained from averaging (14) 
over d> : 



dr 



sgn(g) idL\ =Q 

mcRfio \ dcf> I <f> 



(54) 



where all the terms are ex cose/), sin</>, or cos <f> sine/) and 
average to zero. Consequently, the non-conservation of 
1$, on average, actually starts at 0(1//3q). 

We consider the same set of limiting cases as in Sect. II. 
In the limit of ultrahigh energy, 0q 3> 1$, the con- 
straint (19) with Eq. (20) becomes: 



^(l-w 2 )//3 =w e ± 



(55) 



^J(uje) 2 + (1 - Lu 2 )[n 2 + 2 W sin0 - (1 - w 2 ) sin 2 9} , 

retaining all powers of n. Figs. @ and ffl show the allowed 
regions for several combinations of E and 9q, where the 
allowed regions are defined by the requirement that lg 
in (54) be real. Note that the allowed regions of Figs, g 
and M axe no longer functions of 9 alone, unlike the case 
of eqs. (24,25). The axial symmetry about the dipole is 
lost, but for moderate values of a), the asymmetry is not 
extreme. For large tilts (9 — ► n), the energy n must be 
somewhat higher than in the zero-tilt case to make the 
entire sphere allowed. If 4> averaging is applied to the 
constraint (19), the axial symmetry about the dipole is 
restored. 



FIG. 8. Allowed regions of real lg in the ultrahigh energy 
limit, for r\ = 0.5. Surface coordinates are colatitude 8 and 
longitude <jj m degrees, (a) 9o = 45°. (b) 9o = 135°. 



FIG. 9. Allowed regions of real 1$ in the ultrahigh energy 
limit, for r/ — 1.0. Surface coordinates are colatitude 9 and 
longitude (f> in degrees, (a) #o = 45°. (b) 9q = 135°. 

The limit where 1$ cannot be neglected, but r\ <C 1$, 
is the ultrahigh P^ case; lg is then zero except within the 
narrow allowed regions. If the tilt is zero, eqs. (26-29), 
the particle orbits are almost infinitely narrow circles at 
constant 0. The motion is O(l//3o) relative to the (f> mo- 
tion. With non-zero tilt and 1$ approximately conserved, 
the quasi-circular orbits occur at the same angle as the 
no tilt case, sin0 = y/l$, for 1$ > 0, with O(l//?o) and 
O (Q sin 0o /A) ) corrections . 

The same argument holds for the small angle case, 
which is related to the ultrahigh P$ limit in the limit 
l,j, — > 0, with rj still negligible. The approximate con- 
servation of 1$ still holds in the limit of large /3q, with 
additional corrections of relative order d>sin0o. 

The magnetic flux enclosed by an orbit is, to leading 
order in 1//3q, just as in Sect. II, eqs. (32,36,37). The 
averaged rotational corrections are of order 0{u 2 ). 



of the sphere are forbidden for large tilt and moderate 
rj. The quantum energy levels are still approximately the 
harmonic oscillator levels given by Eq. (39), depending 
only weakly on the field. 

The ultrahigh P$ case results in the energy levels: 

E/mc 2 ~ |l - (w 2 /2)(sin 2 O + l# cos 2 O ) + 



(1 - (5<I> 2 /4) sin 2 O + 01% cos 2 O ) x 



(n e + l/2)^ / 3 /(2mci? v /^)} 



1/2 



(56) 



to leading order in 0q and u>, with non-zero tilt. Note that 
the rotational corrections enter at 0(<D 2 ). The corrections 
multiplying /3o can either raise or lower the energy, but 
the first correction (the pure centrifugal effect) always 
lowers the energy. 

The localized pole states are also corrected by rota- 
tional effects. Their levels arc: 



E/mc 2 ~ |l - (d> 2 /4) sin 2 O + (57) 
(1 - (3w 2 /4)sin 2 o )(n fl + l/2)h/3 /(2mcR)\ 



1/2 



and reproduce Eq. (35) if Q — ► 0. For both this and 
the ultrahigh P^ cases, the rotational corrections begin 
at O(oj 2 ) because of the angular averaging. However, 
corrections of O(u>/f3o) may enter. 

For moderate rotational velocities uj = VtR/c jj 0.1, 
the azimuthally-averaged corrections are 0(uj 2 ). Correc- 
tions not explicitly treated are: O(l//3o), 0(oj/(3q), and 
O(uism9o/(3o) (Fig. no). The multiplicative rotational 
corrections to the magnetic field in the energy arise as 
an effective electric field induced, as seen by an inertial 
observer, by the rotation. Purely rotational (centrifugal) 
corrections also enter. 



5.3 Semiclassical quantization 



In the present case, the quantization of Ia, in the 
Wilson- Sommerfeld condition (33) is only approximate, 
but valid through O(l/f3 ). 

The ultrahigh energy case is not qualitatively different 
if 0q 7^ 0. As seen from the classical orbits, some regions 



FIG. 10. Phase space for the ultrahigh P<p case, with 
negligible rj. Effects of O(l//3o) are greatly exaggerated, (a) 
Phase space for <f) (degrees), showing small variation of I $ for 
the case 8q 7^ 0. The allowed region (phase space trajectory) 
itself has finite O{l/j3o) thickness (not shown), (b) Phase 
space for 9 (degrees) showing small variation of 6 and small 
values of lg . 
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6. CONCLUSION 



We have examined only the simplified scenario of 
charged particles on a rotating spherical surface with an 
intense dipole magnetic field. The resulting energy lev- 
els and magnetic flux quantization, for very high field 
strength, are roughly similar to the relativistic Landau 
and quantum Hall systems. The effect of realistic rota- 
tion velocities (VtR/c ~ 0.1 or less) is small to moderate 
in the two-dimensional case. For these small values of 
u), the case of a tilted axis is not qualitatively different, 
although non-zero tilt creates a non-trivial geometry and 
complicates the dynamics. 

There are a number of related issues in this idealized 
situation not treated here. A treatment of the three- 
dimensional orbitals requires inclusion of the radial mo- 
tion and radial dependence of the vector potential A 
and has been considered for particles confined within the 
crust in a separate work |5| . A fully quantum treatment 
of the orbitals necessitates the Dirac equation and the in- 
clusion of spin, although the basic features of the Dirac 
orbitals are expected to be outlined by the W-S method. 

Application of our results to a neutron star demands 
further realism. In particular, the effect of the lattice of 
nuclei on non-localized electrons and nuclei must be in- 
cluded P,|7|,^2[ . The lattice structure of a neutron star 
crust must be disordered in its surface layers, although 
the fermion temperature is small compared to the Fermi 
energy M. These lattice effects modify the magneto- 
sphere charge cloud, the conduction bands, and the flux 
quantization. The lattice must be included to give a full 
picture of the crust's charge and current distribution. 
The magnetic field itself is also modified if the currents 
contribute significantly to the field. 

Finally, observational questions remain. These include 
verification, if possible, of flux quantization and determi- 
nation of what roles it and the special conduction band 
structure play in the formation and evolution of a neu- 
tron star's magnetic field and crust. 
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